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We show that, in stellar core plasmas, the one-body momentum distribution function is strongly 
dependent, at least in the high velocity regime, on the microscopic dynamics of ion elastic collisions 
and therefore on the effective collisional cross sections, if a random force field is present. We take 
into account two cross sections describing ion-dipole and ion-ion screened interactions. Furthermore 
we introduce a third unusual cross section, to link statistical distributions and a quantum effect 
originated by the energy-momentum uncertainty owing to many-body collisions, and propose a 
possible physical interpretation in terms of a tidal-like force. We show that each collisional cross 
section gives rise to a slight peculiar correction on the Maxwellian momentum distribution function 
in a well defined velocity interval. We also find a possible link between microscopical dynamics of 
ions and statistical mechanics interpreting our results in the framework of non-extensive statistical 
mechanics. 
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I. INTRODUCTION 



In a dense neutral plasma like, for instance, the as- 
trophysical plasma of a stellar core with mean potential 
energy of about the same order of magnitude as thermal 
energy P, 13 j a random microscopical electric field (usu- 
ally called electric microfield) arises at any spatial point. 
Its origin does not rely on correlations between particles 
in the plasma, and indeed it is present also in ideal plas- 
mas of statistically independent particles; rather, it is due 
to local thermal fluctuations in the position of ions 
Commonly, the microfield strength is not negligible, be- 
ing of the same order of magnitude of the Coulomb field 
of a unit charge at the characteristic Wigner-Seitz radius. 
As Romanovsky and Ebeling |^] pointed out, the dynamic 
enhancement on nuclear fusion rates due to electric ran- 
dom fields is large only in very dense stars like white 
dwarfs; on the contrary, its importance inside the Sun's 
core is presently believed to be limited. 

We show that a random field of generic nature (random 
electric or magnetic microfields belong, among others, to 
this category) may play a crucial role, as it influences the 
upper tail of the one-body stationary momentum distri- 
bution function of ions in a dense neutral stellar plasma, 
leading to slight deviations from a pure Maxwellian dis- 
tribution. Furthermore, different elastic collisional cross 
sections among interacting ions could influence signifi- 
cantly the tail of the distribution, and each one provides 
corrections in a characteristic velocity range only. 



'Electronic address: fabrizio.ferroOpolito 
t Electronic address: piero.quarati@polito. 



We want to stress here that, in astrophysical plas- 
mas, many different collisional processes could be ac- 
tive at the same time, provided that we consider dif- 
ferent velocity intervals. Besides, for instance, the usual 
pure Coulomb interaction (described by the well-known 
Rutherford cross section formula), many different screen- 
ing potentials have been proposed since long time, in 
order to provide effective models for several astrophys- 
ical conditions |^. Therefore, screening and many-body 
effects, whose importance relies on the fact that they 
strongly enhance thermonuclear reaction rates, are also 
important, lying on a kinetic framework, as they modify 
the collisional cross section between ions @. 

In this paper we analytically derive the one-body dis- 
tribution function of momentum starting from a kinetic 
equation in which we set three cross sections of inter- 
est in dense and weakly non-ideal plasmas. Our calcula- 
tions are based on the existence of a random force field 
F, which can be justified either by the theory of elec- 
tric microfields (as outlined in Refs. 3, 4]), or the theory 
of dissipative random forces in the approach of Langevin 
equation for Brownian motion (see, for example Ref. [2]), 
and can be originated from density fluctuations in the 
plasma. We are urged in this line of research also by 
the Einstein's criticism to the Boltzmann probability re- 
lation, based on the argument that statistical mechanics 
might only be justified in terms of classical or quantum 
microscopic dynamics Q]. 

In the first part of our paper (Sections Ull IIIII and II V|) 
we use a kinetic approach to describe the motion of parti- 
cles submitted to a generic random force field with finite 
relaxation time (i.e. not (5-correlated) . The stationary 
solution can be expressed in terms of collisional cross 
sections and collisional frequencies; therefore we may es- 
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tablish a link between the type of particle collisions and 
the form of stationary distribution functions that can dif- 
fer from the equilibrium Maxwellian function. Then, we 
define the parameter q characterizing the deformation 
factor of our distribution and we calculate it in terms 
of known physical quantities. Interpreting this deforma- 
tion on the physical ground of non-extensive statistical 
mechanics, and Tsallis statistics as a special case 
the parameter q can be expressed in terms of dynami- 
cal quantities like cross section, ion-correlation parame- 
ter and plasma parameter. A dynamical realization of 
non-extensive statistics and, in general, of superstatis- 
tics has been completed by Beck [l^l by using stochastic 
differential equations with spatio-temporally fluctuating 
parameters. 

In the second part (Section 0), we consider the quan- 
tum energy-momentum uncertainty due to contempora- 
neous interaction among many particles of plasma. We 
use the Kadanoff-Baym ansatz ,11.1 a nd follow the ap- 
proach of Galitskii and Yakimets [l^l that leads to the 
appearance of power-law tail in the momentum distri- 
bution function. This function shows an enhanced tail 
at high momentum depending on the collision frequency 
and the coUisional cross section. As a consequence, we 
show a link between non-extensive statistical mechan- 
ics and dynamics, by evidentiating the type of micro- 
scopic elastic collisions acting among particles and how 
the collision frequency is related to momentum distribu- 
tions. Then we investigate the cross section that repro- 
duces the non-extensive distribution function, limiting 
ourselves, for simplicity, to the case q > 1, and we show 
that the corresponding interaction looks like a tidal-like 
force superimposed on the two-body attraction, giving a 
collision cross section a{£p) ~ where Sp is the rela- 
tive kinetic energy. The quantum approach that we fol- 
low along our discussion relies on equilibrium conditions 
of the system; therefore our final result is a real equilib- 
rium (not metastable or stationary) distribution function 
that differs from the equilibrium Maxwell-Boltzmann dis- 
tribution. The existence of such an unusual equilibrium 
function is due to the quantum uncertainty as it is shortly 
described in Section Ivl 



II. KINETIC EQUATION UNDER A 
GENERALIZED RANDOM FORCE 

A kinetic equation describing electrons in a plasma 
under an external electric field, in which collisions be- 
tween electrons and neutral atoms are present, was de- 
rived by Chapman and Cowling ITp|. bv Spitzer and 
by Golant, Zilinskij and Sacharov |T3- They express the 
actual one-body velocity distribution function as a formal 
series, 

/(V) = / + /i + /2 + . . . , 

where / = f{v) is the isotropic component (with v = |v|), 
while /i,/2 and so on describe next orders of anisotropy 



induced by the external field. In addition, f{v) may be 
only a slight perturbation of the Maxwellian distribution 
function. 

Here, we adopt their equation in order to derive the 
momentum stationary distribution of ions, but we replace 
the external electric field with a generalized random force 
F and we focus on the isotropic function / only. The 
elastic coUisional cross sections that we are considering 
describe the interaction among ions and among ions and 
electric dipoles of polarized neutral compounds of the 
Wigner-Seitz spheres. All these cross sections will be 
discussed in the next sections. 

Thus, considering the plasma component consisting of 
ions of mass m, the kinetic equation reads 0| 



^2 F2 



d/ 



3 /x^J/^ d?; 



ksTd/ 



A* 
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where v is the modulus of the relative velocity between 
two ions, fi is their reduced mass, i^iv) is the coUisional 
frequency, keT is the thermal energy and k, which is 
defined as 

. = 24, 

is the energy transfer coefficient (or an average value of 
it). 

Let us discuss briefly the origin of the double sign in 
the term containing the random force F in Eq. ^ The 
quantity 



Df 



2 F2 
±-- 



3 fi^v'^ 



is the perturbation on the diffusion coefficient of the sys- 
tem due to the force F, and the corresponding particle 
fiux is given by 

Jf = tDfti— , 
av 

n being the ion particle density of the plasma. If F were 
an electric microfield, i.e. F = eE (with e equal to the 
electric charge of one ion), the corresponding sign would 
be positive, thus enhancing the actual diffusivity of the 
system, while in the opposite case, the total diffusivity 
would drop. Therefore, we introduce the double sign 
since we want to deal with the most general situation, 
in which either sub-diffusivity or super-diffusivity may 
be significant. 

The analytical solution of Eq. ^ reads 



f{v) = /(O) • exp 



(2) 



where the constant /(O) should be calculated through the 
normalization condition 



+00 



Av'v'" f{v') — 1 
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Let us define now a characteristic strength of the gen- 
erahzed random field F as 



Fc = vy/ K^keT 



Then, if the condition 



F^ < 



c > 



(3) 



holds, i.e. if the random force is neghgible, Eq. |21 gives 
the Maxwellian distribution function at temperature T. 
The central point is that, in this case, the Maxwellian 
function is the solution of the kinetic equation ^ regard- 
less of any assumption about the coUisional frequency v 
of the plasma. 

On the contrary, if the condition in Eq.|31fails, the form 
of the solution f{v) is determined by the explicit depen- 
dence of the coUisional frequency v on relative velocity v. 
The i'{v) frequency is itself a function of the coUisional 
cross section (t{v), as 



= nvcr(v) . 



(4) 



Thus, in this case, Eq.|21 leads to a Maxwellian pro- 
vided that we choose the cr{v) = aov~^ cross section (ao 
being a suitable constant), and that we renormalize the 
temperature of the plasma in the following fashion, 



^bTcS = keT ± 



F' 



3 Kfin^aQ 



1 2 ' 



(5) 



where Tos is an effective temperature which will be of 
central importance in our subsequent discussion. 



III. ANALYTICAL SOLUTION OF THE 
KINETIC EQUATION 

In the following, we discuss the effect of three different 
cross sections, ctoj fi and tT2, whose explicit functional 
dependence on relative velocity, together with that of 
their coUisional frequencies, is respectively. 





= a^v ^ 






= nao , 


(6) 
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= na-\V , 


(7) 
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(8) 



where ao, ai and a2 are dimensional constants. We shall 
discuss in Section llVl the physical meaning of the previous 
cross sections in dense astrophysical plasmas. 

We state the hypothesis of absence of interference be- 
tween the three collision types, namely we assume that 
the total coUisional frequency v could be cast in the fol- 
lowing approximate fashion. 



because different types of collision act significantly only 
in separate velocity intervals as evident from the func- 
tional dependencies reported in Eqs.|ni[7|and|Sl 
Let us express now the solution of Eq. |2| as 



fiv) = /(O) • exph/(«)] , 
where we have defined the integral function 

fiv'dv' 



(9) 



liv) 



keTil-j;^ 



3 K. fin^a^ l+cii;'^+C2f 

v'dv' 



1 



(10) 



l + Cit)'^+C2U'* 

\2 



with ci = (ai/ao)^, C2 = (a2/ao)'' and r = Tes/T - 1, 
according to Eqs.[Slini[71and|Sl 

From Eq. 1101 we immediately obtain 



I{v) 



2kBT " 2kBT 



Th(v) 



where 



hiv) 



du- 



1 



/O C2U^+CiU + T+l 

Let us define now the following parameter. 



(11) 



(12) 



K = -■ 



-1 

4C2 



1 



whose sign is physically relevant as we are about to show. 

If _ftr < 0, Eq. 1121 gives, apart from an imimportant 
numerical term. 



h = 



1 



2C2V^ + Cl - 2y/\K\^ \ 

In ^^^^= -I- const 



2C2V^ 



Cl 



C2 



which in turn, through Eqs. 1111 and El yields our first 
result (for K < 0) 



f{v) oc exp 




In this case, the cross section di dominates and the 
generalized random force field F does not play any role in 
the region of interest for fusion reaction rate calculations 
in astrophysical plasmas, as the perturbation from the 
Maxwellian function vanishes in the limit v —^ +00 (nev- 
ertheless, it can be of interest in studies of some atomic 
processes like radiative recombination, whose cross sec- 
tion increases as v goes to zero and that, therefore, has 
rates sensibly modified by slight corrections at low veloc- 

ity). 
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As far as astrophysical plasmas are concerned, a 
more interesting physical situation occurs when K > 0, 
namely, if the condition 



Aalal ' 



(13) 



holds. Condition 1131 is fulfilled providing the force field 
F is 



/3 , ^( a\~ Aalal 
F>nJ-.^k^Ti 



in case of super-diffusivity, or instead 



F <n^ -K/ikeT ^ 

when considering sub-diffusivity. 
From Eq. E|we obtain 



(14) 



(15) 



h{v) = 



1 

K 



du 



ci 



K 2VK^ 



— arctan - 



C2 2 

arctan | ^ / — v 
K 



Cl 



ci 



(16) 



Starting from Eos. 1111 1121 and [TCI we can express I{v) 
as a formal series of powers of w^, 



I{v) 



where 



and 



2kBTeff 



2kBTc 



cff 



+ 7 



2kBTeff 



9 F2 2 

5 = ±-- "1 



3 2 2 4 ' 



7 = ±7T 



8 F^keT , 



9 K^^n? Q!q 



+ 



both being \5\, I7I < 1. 

Therefore, the final form of the one-body distribution 
function under generalized random fields reads 



fiv) cx exp 



2kBreff 



X exp 



/it) 



exp 



-7 



that corresponds to an energy probability factor 



f{ep) cx exp 



X exp 



-S 



keTcff 



exp 



(17) 



where £p = p^/(2/i) is the centre-of-mass kinetic energy, 
given the linear momentum p = /xv. 

It is noteworthy that our result in Eq. El ma-y be re- 
lated at least for small deformations to the non-extensive 
distribution function at the same effective temperature 
Tos 9] 



/(ep) « 



(18) 



where q is called the entropic parameter. As can be 
straightforwardly shown after some manipulations, in the 
low deformation hmit {q — l)£p/(kBroff) 0, Ea.lT^re- 
duces to Eq. El provided that d = {1 — q)/2. Thus, this 
condition establish a link between the entropic param- 
eter q and our parameter d that comes from microfield 
strength and cross sections. We point out that, in the 
same limit, also other distributions of generalized statis- 
tics have the same behaviour, as explained in Ref. ^2^. 

We recall that in recent past we have shown that, if the 
generalized random force is due to an electric microfield 
distribution, the parameter S can be related to the non- 
extensive (Tsallis) entropic parameter q and the following 
analytical expression can be derived. 



2„4 



6 = = Ur^a 

2 

where F is the plasma parameter and a is a dimension- 
less parameter accounting for ion correlations in the ion- 
sphere model (0.4 < a < 1) 

From Eq. El it follows that there exist three dif- 
ferent intervals of relative velocity in which the evalu- 
ated corrections due to the random field are sufficiently 
large to be noteworthy. First of all, we observe that if 
Ep ~ keTefi, the dominant factor is exp[— ep/(kBrcff )], 
namely the Maxwellian factor characterized by the cross 
section cfq. The exponential factor with the S parameter, 
corresponding to the coUisional cross section ai , becomes 
not negligible with respect to the Maxwellian only when 
£p ^ kBTcff/l^l; it is also often called the Druyvenstein 
factor. Finally, the third term exp[— 7(ep/kBToff )'^] arises 
when Ep ^ |5/7|kBToff; as we shall sketch in Section^ it 
can be related with quantum energy-momentum uncer- 
tainty effects in dense astrophysical plasmas. 

In conclusion of this section, let us summarize that if 
the random force field is absent or negligible, in spite 
of the presence of whatever kind of coUisional cross sec- 
tions, all stationary states which are solutions of the ki- 
netic equation have an analytical expression that coin- 
cides with the equilibrium Maxwellian distribution func- 
tion. Therefore, the non-extensive statistical description, 
at least in a classical framework, requires as a first con- 
dition that particles be subjected to a sufficiently strong 
random force field and, as second condition, that a con- 
stant coUisional cross section (or depending on higher 
positive powers of velocity) should act among the parti- 
cles of the system. 
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IV. COLLISIONAL CROSS SECTIONS IN 
ASTROPHYSICAL PLASMAS 



We want to discuss the physical meaning of the coUi- 
sional processes related to ctq, cti and (T2, defined in Sec- 
tion lllll and inserted into the kinetic equation in order to 
derive the one-body distribution function. 

The cross section (7o{v), defined in Eq. El is the most 
important one as it originates the well-known Maxwellian 
distribution function even in presence of a generalized 
random field. That the solution of the kinetic equation 
(Eq. nj, at first order, be the Maxwellian function is our 
first unavoidable requirement, because the actual kinetic 
solution for F = is indeed the Maxwellian, and we are 
dealing only with slight corrections. 

Following Ref. ^3 we can state that, starting from an 
interaction force that depends on distance as r~^, the 
corresponding cross section is cr{v) (x v~^^^'^~^\ There- 
fore, in the case of the cross section (To{v) oc v^^, the 
underlying force goes like r~^, while the potential en- 
ergy is proportional to r~^, and we can interpret it as 
the cross section due to the interaction between an ionic 
charge and an electric dipole induced by the ion on the 
neutral system of charges composing a Debye sphere jlSj . 

On the contrary, if we considered a pure Coulomb in- 
teraction due to a force Fc oc (with s = 2), it would 
give a coUisional cross section proportional to v~'^; how- 
ever, this case is not physically suitable in presence of 
an intensive random force field because of induced di- 
vergences in the distribution function at small velocities. 
Krook and Wu showed in the past that coUisional 
cross sections going like and always produce, 
after a sufficient long time, a Maxwellian distribution; 
however their system is not subjected to a random force 
field. 

The cross section o'i(f), defined in Eq. [7\ has been 
introduced by Ichimaru 19] in the framework of an ion- 
sphere model for non-ideal and weakly-coupled plasmas 
with a r parameter of order unity, and with a small num- 
ber of ions inside the Debye sphere. In these physical 
conditions the coUisional cross section, directly derived 
from the pure Coulomb one, is constant according to the 
approximations of the model and it can be cast in the 
simple form 

(Ji{v) w 27r(Q;a)^ , 

where a is the correlation factor already introduced in 
Section IlIII and a is the interparticle distance. The cor- 
rection due to fJi on the probability function of energy 
is of order exp[— (5ep/(kBToff )^], and shows the same be- 
haviour of the so called non-extensive corrective factor 
(see Ref. (§3). 

The cross section a2{v) will be described in the next 
section. 



V. EFFECTS OF QUANTUM 
ENERGY-MOMENTUM UNCERTAINTY ON 
THE EQUILIBRIUM DISTRIBUTION 
FUNCTION 

Here we introduce simple arguments to synthetically 
explain the meaning and justify the use of the cross sec- 
tion a2(v) defined in Eq.|Sland, at the same time, to show 
a possible link between quantum energy-momentum un- 
certainty and non-extensivity. 

Quantum energy-momentum uncertainty in weakly 
non-ideal dense stellar plasmas influences thermonuclear 
rates. In fact, in classical physics, energy e and mo- 
mentum p (or Ep = p^/2yU, with n reduced mass) are 
linked by the dispersion relation S{e,ep) = S{e — Sp). 
Nevertheless, if the particles cannot be considered free, 
e and Sp are independent variables. For instance, an ion 
tunnelling the Coulomb barrier during a thermonuclear 
fusion reaction can interact simultaneously with many 
other particles. In this case, the dispersion relation is 
given by the function 6j(e,Sp) defined using the ansatz 
of Kadanoff and Baym jlj. Under equilibrium condi- 
tions, and this time without any random force field, the 
energy and momentum generalized distribution function 
can be written as 

/(e,£p) = o^(e,ep) , 

TT 



with 



ImE«(£,£p) 



{e-Ep- Rel]«)2 -I- (Iml]^)2 ' 



where n(e) is the particle number distribution and 
T,^{e,ep) is the mass operator for the one-particle re- 
tarded Green function. 

Galitskii and Yakimets [13] (see also Refs. jH 12^ ) 
derived that the quantum energy-momentum indetermi- 
nacy and a non-zero value of ImS^ lead to the non- 
exponential tail of the energy probability factor f{sp). 

We limit ourselves to the case of a dispersion relation 
of Lorentz type, 

5 is e \ -- ^(^'^P) 



with 



7(e,ep) = hv'l^u(e,£p) = hnatisp) 



2e 



1/2 



where v'^oiA^iSp) is the total collision frequency and 
(Jt{£p) is the coUisional cross section. 

Let us make the example of a pure Coulomb interac- 
tion. We have that 



7(£,ep) 



6 



and the momentum distribution becomes 



with 



and 



7r3/2 (kBT)3/' 



fMB{£p) = exp 



[fMBiep) + fQ{Sp)] ,(19) 



fQiSp) = f C(kBT)3/2 1 , 

where C is a constant depending on the density n. 

At high momenta the last term in Eg. 1191 can be many 
orders of magnitude greater than the first one and rep- 
resents an enhancement of the tail, with important con- 
sequences in the calculations of nuclear fusion rates. 

We want to verify if, by using a certain elastic col- 
lision cross section, we could obtain from the quantum 
uncertainty effect the non-extensive Tsallis distribution, 
limiting ourselves, for simplicity, to the case of entropic 
parameter g > 1 

Following (and adapting to our present needs energy 
fluctuations instead of temperature fluctuations) the ap- 
proach outlined by Beck and Cohen |2^, we may state 
that any non-Maxwellian energy probability function 
should be cast in the form of a Laplace transformation 
of the function (5 f(-E, £p) which describes the non ideality 
of the system |2J|, 



r+oo 

f{ep) = I dEexp 



E 



5f{E,ep) 



The function 5f{E,ep) must be assumed to be a 
Gamma- (or x^-) function, in order to obtain that f{sp) 
be a non-extensive (Tsallis) distribution 10, 23]. 

Let us compare this integral with the integral of 
Eq. El that can be written expHcitly as /(£p) = 
/ de exp(— e/kBT)5-y(e, Sp). Quantum uncertainty and 
non-extensivity are two different and distinct causes 
of deformation of the Maxwell-Boltzmann distribution. 
Nevertheless they give the same effect if the microscopic 
interaction among the particles (i.e. the coUisional cross 
section) is of a particular nature as we discuss below. 

Let us pose attention to the physical property of inter- 
est is the 'width' Dj of the Sf{E,£p) distribution; it can 
be shown that for the non-extensive distribution func- 
tion, the width is Dj^e ^ £p, while for the quantum 
uncertainty the Dq ^ a'^{ep)£p relation holds, where 
a{ep) is the collisional cross section. If we now impose 
that super-extensivity and quantum uncertainty give the 
same physical effect on distribution functions, we should 
require that cr(£p) cx or, in terms of relative velocity, 
cr[v) = 172(1') oc V. Thus, the cross section tT2 that we 
have used in Section IIIII is strongly related with both 



quantum and non-extensive statistical effects. The non- 
extensive and the Galitskii-Yakimets distributions result 
given by the same expression. 

Let us recover the behaviour of the interaction force 
responsible of the cross section cr(ep) ~ We can 

write its dependence on the relative coordinate r of the 
two interacting particles as 



Fir) = /o 



Rr 



where /o is a dimensional constant, Rq is a characteristic 
distance of the two-body center of mass with respect to 
a given origin and s is a negative or positive integer. 
Defining the collisional cross section as 



a ~ nd 



with 



JO] — rj- 



in order to have the requested behaviour of cr(£p) oc 
we must set s = —3. Let us recall that the case s — —3 is, 
from the point of view of the orbit differential equation of 
motion, one of the integrable cases, with solutions given 
in terms of elliptical functions . 

Therefore the interacting force responsible of the col- 
lisions that lead to cr(ep) ^ y/s^ reads 



r<Ro 
r> Rq 



where the cut-off is needed in order to avoid divergences 
of the potential energy. 

We may argue that the force Fq (r) can be understood 
as a tidal-like force 26] if we assume that an attractive 
central force of intensity /qq, centered at a distance Rq 
from the center of mass of the two interacting particles 
separated by a distance r, is superimposed. The tidal- 
like force acts globally over all the particles of the sys- 
tem. This is the dynamical requirement to recover the 
non-extensive (Tsallis) distribution in the framework of 
quantum energy-momentum uncertainty. It is worth and 
peculiar to remark that by applying the virial theorem to 
this case we obtain a negative kinetic energy, which is, in 
fact, understandable and admissible by the uncertainty 
principle. 

We derive the analytical expression of q by equalling 
the complete expressions of Dme and Dq. We obtain 



9=1 + 



{hcfnH^iY Rl 



where n is the plasma density and is of order of 
unity p^ ]. 

The correction to the unity can be thought due to the 
many-body effect over the two-body interacting system. 
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As an example, let us make the following numerical ap- 
proximate evaluation of /q^ : if the correction is of order 
of ten per cent, the density n 10^^'* fm^'', and Rq « 
10^ fm, we obtain, for a proton plasma (/ic^ w 460 MeV), 

/q„ « 10-12 MeV/fm. 

Before concluding this section we remark that the 
non-extensive distribution usually describes metastable 
states or stationary states of non-equilibrium systems. 
On the contrary, in this case, quantum uncertainty with 
coUisional cross section cr(£p) ~ gives a distribu- 

tion function which belongs to an equilibrium state, al- 
though different from the Maxwell-Boltzmann distribu- 
tion. Other generalized distributions have been recently 
proposed 27] . For situations with small deformation our 
argument are valid also for these distributions. 

VI. CONCLUSIONS 

We have set a kinetic equation suitable to describe the 
stationary states of a weakly non-ideal plasma of a stel- 
lar core subject to generalized random forces. Provided 
that a random force satisfying condition ^] for super- 
extensivity or 1151 for sub-extensivity is present, the mo- 
mentum distribution function can be cast in the sim- 
ple fashion of Eq. 1171 in which, besides the well-known 
Maxwellian factor, other terms also contribute. 

The momentum distribution function is formally iden- 
tical to the non-extensive distribution (when q > I) ex- 
panded in powers of (1 — g) for small deformation. An 
analytical expression of q, the entropic parameter, can 
be derived in terms of the elastic collision cross sections 
acting among the particles of the system. 



The main point is that each correction factor is due 
to a particular coUisional process between ions, and that 
each of them contributes in a well defined interval of rela- 
tive velocity, as shown at the end of Section llTTl All these 
corrections are small, nevertheless they are not negligi- 
ble at high energy, i.e. in the region of ion spectrum of 
predominant interest for calculations of nuclear reaction 
rates in astrophysical plasmas. 

We have stressed that in physical conditions as, for 
example, stars with F > 1, many coUisional processes 
may be active, maybe at the same time, and that each 
one of them is described by a cross section with a de- 
pendence over velocity stronger (proportional to v~^, 
or even cx v) than the simple Coulomb scattering (pro- 
portional to v~'^). This fact is intimately related to sta- 
tistical many-body effects and represents a link between 
dynamics (the type of two-body elastic collisions) and 
statistical mechanics (the momentum distribution func- 
tion of the stationary states involved). 

Finally, in the framework of a quantum many-body de- 
scription of the equilibrium state, considering the energy- 
momentum uncertainty due to the non-commutativity of 
position and momentum operators, we have found that 
if the coUisional cross section cr(ep) behaves like ^/e^, 
the distribution function coincides with the non-extensive 
(Tsallis) distribution function with q > 1. The requested 
behaviour of the cross section (j{sp) is due to an inter- 
action similar to a tidal-like force. Therefore the anal- 
ogy between quantum uncertainty effect on the distribu- 
tion and non-extensive effect is achieved provided that 
an overall attractive interaction is superimposed to the 
two-body interaction. This represents again a possible 
link between dynamics and statistical mechanics. 
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